Polarization of the prompt 7-ray emission from the 
7-ray burst of 6 December 2002 
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Observations of the afterglows of 7-ray bursts (GRBs) have revealed that they 
lie at cosmological distances, and so correspond to the release of an enormous 
amount of energy— The nature of the central engine that powers these events 
and the prompt 7-ray emission mechanism itself remain enigmatic because, once 
a relativistic fireball is created, the physics of the afterglow is insensitive to 

X' 

the nature of the progenitor. Here we report the discovery of linear polariza- 
tion in the prompt 7-ray emission from GRB021206, which indicates that it is 
synchrotron emission from relativistic electrons in a strong magnetic field. The 
polarization is at the theoretical maximum, which requires a uniform, large-scale 
magnetic field over the 7-ray emission region. A large-scale magnetic field con- 
strains possible progenitors to those either having or producing organized fields. 
We suggest that the large magnetic energy densities in the progenitor environ- 
ment (comparable to the kinetic energy densities of the fireball), combined with 
the large-scale structure of the field, indicate that magnetic fields drive the GRB 
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explosion. 

We used the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI)— to 
make these7-ray observations of GRB021206. RHESSI has an array of nine large-volume 
(300 cm 3 ) coaxial germanium detectors with high spectral resolution, designed to study solar 
X-ray and7-ray emission (3keV-17MeV). RHESSI has high angular resolution (2") in the 
~1° field of view of its optics; however, the focal plane detectors are unshielded, open to the 
whole sky. Thus, while the chances are small that RHESSI will see a GRB in its imaging 
field of view, it measures them frequently in the focal plane detectors themselves. These 
observations provide highresolution spectra, individual photon times and energies, as well 
as the potential for polarization measurements. RHESSI is not optimized to act as a 7-ray 
polarimeter, but several aspects of its design make it the most sensitive instrument to date 
for measuring astrophysical 7-ray polarization. 

In the soft 7-ray range of ~0. 15-2.0 MeV, the dominant photon interaction in the RHESSI 
detectors is Compton scattering. A small fraction of incident photons will undergo a single 
scatter in one detector before being scattered and/or photoabsorbed in a second separate 
detector, which are the events sensitive to the incident 7-ray polarization. Linearly polar- 
ized 7-rays preferentially scatter in directions perpendicular to their polarization vector. In 
RHESSI this scattering property can be used to measure the intrinsic polarization of as- 
trophysical sources. The sensitivity of an instrument to polarization is determined by its 
effective area to scatter events, and the average value of the polarimetric modulation factor 

n(6, E), which is the maximum variation in azimuthal scattering probability for polarized 
photons. This factor is given by fi = (dcr± — do~\\)/(do~± + do~\\), where do~±,do~\\ are the Klein- 
Nishina differential cross-sections for Compton scattering perpendicular and parallel to the 
polarization direction, respectively, which is a function of the incident photon energy E~, 
and the Compton scatter angle 6 between the incident photon direction and the scattered 
photon direction. For a source of count rate S and fractional polarization IL, the expected 
azimuthal scatter angle distribution is dS/dcj) = (S/2tt)[1 — /x m II s cos(2(</> — 77))], where (f> is 
the azimuthal scatter angle, rj is the direction of the polarization vector, and \i m is the aver- 
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age value of the polarimetric modulation factor for the instrument. Although RHESSI has 
a small effective area (~20cm 2 ) for events that scatter between detectors, it has a relatively 
large modulation factor in the 0.15-2.0 MeV range, /i m ~ 0.2, as determined by Monte Carlo 
simulations described below. 

In comparison with other 7-ray instruments (COMPTEL, BATSE) that have attempted 
to measure polarization in the past>^, RHESSI has the major advantage of quickly rotating 
around its focal axis (centred on the Sun) with a 4s period. Rotation averages out the 
effects of asymmetries in the detectors and passive materials that could be mistaken for a 
modulation. Because polarimetric modulations repeat every 180°, any source lasting more 
than half a rotation (2 s) will be relatively insensitive to the systematic uncertainties that 
typically plague polarization measurements. Finally, although the RHESSI detectors have 
no positioning information themselves, they are relatively loosely grouped on the spacecraft, 
allowing the azimuthal angle for a given scatter to be determined to within A0 = 13° r.m.s. 
This angular uncertainty will decrease potential modulations by a factor of 0.95, which is 
included in our calculated modulation factor. 

Prompt 7-ray emission from GRB021206 was detected with RHESSI on 6 December 
2002 at 22:49 UT (Fig. [TJ. This GRB was also observed^ with the Interplanetary Network 
(IPN), which reported a 25-100 keV fluence of 1.6 x 10~ 4 ergcm -2 , and a peak flux of 2.9 x 
10~ 5 ergcm -2 s _1 , making this an extremely bright GRB. The IPN localized 8 GRB021206 to 
a 57 square- arcminute error box located 18° from the Sun. For RHESSI we analysed photons 
in the energy range 0.15-2.0 MeV that scattered between two, and only two, detectors for the 
5 s integration period shown in Fig. [I] Scattered photons constitute roughly 10% of the total 
0.15-2.0 MeV light-curve events. Counts were binned by the centre-to-centre azimuthal angle 
between the two detectors around the RHESSI roll axis, corrected for the rotation of the 
spacecraft at the time of the photon event. This azimuthal distribution is plotted in Fig. El 
The top panel shows the raw data, as well as the expected variation for an unpolarized GRB 
due only to the light-curve variability. The bottom panel shows the residual of the data 
once this unpolarized distribution is subtracted. The residual shows a large modulation, 
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which we interpret as a linear polarization of IT m = (80 ± 20)%. This observation is the 
first astrophysical polarization measurement at 7-ray energies. Note that the uncertainty on 
this polarization amplitude reflects in large part our uncertainty in the modulation factor. 
The fact that we have measured a polarization somewhere in this range has a significance of 
< 10 -8 (or a confidence > 5.7a). 

Linear polarization is generally considered a clear indication of synchrotron emission. For 
electrons with an energy spectrum characterized by a power-law distribution with spectral 
index p, synchrotron photons are emitted with a linear polarization- of II = (p + l)/(p+7/3). 
For shock acceleration— , typical values of p =2-3 correspond to linear polarizations of 70- 
75%. For unresolved sources, polarizations from many directions generally add together to 
produce net polarizations that are a fraction of this maximum value. While the source is 
unresolved in our observation, if the source electrons are moving with a bulk Doppler factor 
of r, we are only viewing the source over a solid angle f2 7 « 1/T 2 . Thus, for typical values of 
T > 300 that have been implied from GRB afterglow observations, we are effectively resolving 
a source region of solid angle fi 7 ~ 10~ 5 sr. Our measurement of this high polarization 
is consistent with synchrotron origin for the initial GRB from a region of nearly uniform 
magnetic field. For this emission process to be radiatively efficient as is implied by many 
afterglow observation s 11 ' 12 , the magnetic energy densities must be close to equipartition 
(comparable to the kinetic energy densities of the fireball) 1 ?i 14 . 

This polarization from the prompt 7-ray emission is significantly higher than optical 
polarizations of 1-3% typically measured from afterglow s 15 i 16 i 17 ; as well as the optical polar- 
ization of 10% recently reported^ from the afterglow of GRB020405. The afterglow emission 
implies a strong magnetic field behind the shocks^, although the field energy density can 
be well below equipartition 20 . However, the implied magnetic field strengths are too large 
to have been due to a progenitor field being dragged along by the expanding fireball—, 
or compression of the ISM magnetic fields in external shocks 22 . Therefore, the magnetic 
fields responsible for the afterglow synchrotron emission are probably turbulent fields that 
have built up behind the shocks 21 :23 , which is consistent with the relatively small optical 
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afterglow polarizations. This locally generated field invoked for the afterglow has influenced 
many researchers to consider the magnetic field responsible for the prompt burst of 7-rays 
as a turbulent, fireball-induced field as well. In the standard internal shock model—, the 
prompt 7-rays are produced by synchrotron emission of electrons accelerated to relativistic 
energies by shock acceleration, requiring magnetic energy densities near equipartition in the 
progenitor environment. However, late-time turbulent fields have no direct implications on 
whether the fields responsible for the prompt 7-ray emission are predominantly turbulent or 
organize d 21 ^ 25 . Our observation conclusively shows that the engine driving the GRB has a 
strong, large-scale magnetic field. 

Another potential source of polarization would be as follows: if unpolarized 7-rays are 
initially beamed into a small-angle jet, and then scatter at an angle 9 into our line of sight. 
The polarization induced by this Compton scattering— is II = (1 — cos 2 #)/(l + cos 2 9) for 
photon energies below ~0.1MeV, decreasing as the photon energy approaches and exceeds 
the electron rest mass, 0.511 MeV. For our energy band, we estimate a maximum Compton- 
induced polarization of 70% for scatter angles near 90°. However, this process would be 
inefficient, requiring a much larger 7-ray energy budget than the synchrotron case. First, the 
distribution of scattered 7-rays would be nearly isotropic, requiring an energy budget Air /VLj 
larger than if the observed 7-rays originated from a collimated jet of opening solid angle 
Qj. Second, to maintain such a high polarization as we observed, the 7-rays must undergo 
only a single scatter into our direction because secondary scatters will erase the induced 
polarization from the initial scatter. This condition requires that the scatter medium be 
optically thin, r << 1, and that the luminosity of the unseen initial 7-ray beam be larger 
than the observed scattered 7-ray emission by a factor, 1/r. Because Compton-scattering- 
induced polarization requires a total 7-ray luminosity several orders of magnitude larger 
than that implied by synchrotron emission, and elaborate source geometries, the synchrotron 
origin for the polarization is preferred. 

We suggest that our observation is evidence that the magnetic fields are actually powering 
the GRB explosion itself. It has been argued that a "passive" magnetic field-that is, a 
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field dragged from the surface of central object with a magnetic dipole moment, but not 
driving the GRB-could not be strong enough to produce the prompt 7-ray emission without 
an additional locally generated turbulent magnetic field^. If this conclusion holds, our 
observation is consistent with models of a magnetically driven GRB fireball, where the 
driving magnetic field was generated by extracting the rotational energy of an accretion disk 
around a central compact object through differential rotatio n 2 ^ 2 ? , by directly extracting the 
spin energy of a black hole threaded by magnetic field lines- 2 ^, or by extracting the spin 
energy of a highly magnetized neutron star—. Alternatively, our observation of a large-scale 
magnetic field could support models of dynamos in the post-shock flows^ if the shocks can 
be shown to be unstable on large size scales and on timescales comparable to the prompt 
7-ray emission. 
Methods 

The top panel of Fig. El shows the azimuthal scatter distribution of 0.15-2.0 MeV events 
for the 5 s integration period shown in Fig. ^ corrected for the rotation of the spacecraft at 
the time of each photon event. This distribution is the sum of three components: the GRB 
scatter event rate which averages 820 ± 8 counts per bin, the chance coincidence rate which 
averages 374 ± 6 counts per bin, and the background scatter event rate of 49 ± 2 counts 
per bin. These rates were determined and confirmed using event rates before and after the 
burst, combined with studies of the readout times of multiple detectors during scattered 
events, and were verified independently using our Monte Carlo simulations. Although the 
rotation will average out systematic variations in the scatter angle distribution, we still have 
to correct for the complex time profile of the burst itself, which will cause variations for 
an unpolarized source owing to the finite number of potential scatter angles RHESSI can 
measure at any given instant. We modelled this effect by using the 0.15-2.0 MeV total count 
rate in the RHESSI instrument (Fig. as the time dependent flux template for a photon 
transport Monte Carlo simulation, and using the time-averaged GRB photon spectrum as 
measured by RHESSI for our input spectrum. This simulation used the detailed RHESSI 
mass model that has been developed under CERN's GEANT package, allowing us to model 
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the instrument response to a GRB at the IPN 8 sky coordinates for each rotation angle and 
instantaneous flux, assuming an unpolarized source. This distribution is also presented in 
the top panel of Fig. El We are looking for a modulation signal relative to this variation 
induced by the GRB time profile. 

In the bottom panel of Fig. |2]we show the residual of the measured distribution once we 
have subtracted away the simulated response for an unpolarized GRB, showing our absolute 
modulation signal. For an unpolarized source we would expect this distribution to be flat, 
which we can rule out at an extremely high confidence level (x 2 = 83.5, 11 degrees of freedom, 
d.f.). When we fit this with a modulation curve, the fit improves significantly (x 2 = 16.9, 
9 d.f.), with an amplitude of 128 ± 16 counts per bin. Statistically, this is a reasonable fit to 
the data (the probability of x 2 > 16.9 is 5%), but could be improved with even more detailed 
Monte Carlo simulations including timedependent spectral variability. Using the count rates 
given above and the simulated distribution for an unpolarized GRB, we performed further 
numerical simulations to determine the probability that an unpolarized GRB could produce 
a modulation as large as the one we measure due to random Poisson counting statistics. 
We found this probability to be very low, < 10~ 8 , which translates to a confidence that 
we have measured a polarization at a level > 5.7a. Finally, we estimated the modulation 
factor to be fi m = 0.19 ± 0.04, using both a separate photon transport code which fully 
treats polarization in scattering and uses a simplified mass model, as well as analytical 
estimates based on the GEANT simulation with the full RHESSI mass model. Combining 
the modulation amplitude, the total source scatter event rate, and the RHESSI modulation 
factor, we derive a measured polarization Il m = (80 ± 20)%. 

A number of tests were performed to check that themeasured modulation is real. First, we 
verified that the simulated variation induced by the GRB light curve is accurate by comparing 
it to an angular distribution of events that were chance coincidences in two detectors. These 
interactions are nearly simultaneous, but separated by enough time to distinguish them as 
chance coincidences, not real scattered photons. This distribution should exhibit the same 
variations owing to the GRB light curve, but no polarization. When we subtracted the 
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simulated distribution from the chance-coincident distribution, we found no evidence for a 
residual modulation. We have performed a number of independent checks to make sure 
we do not see modulations from other sources as well. We have verified that extended 
RHESSI background observations show no sign of modulations. In addition, we have done 
a preliminary analysis of a strong solar 7-ray flare observed on 23 July 2002, where we see 
some evidence for a modulation, but corresponding to a polarization < 10%. Therefore we 
feel confident that we have characterized the systematic effects in RHESSI to below the 10% 
polarization level. 
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Figure 1: RHESSI light curves (in total measured counts) in three energy bins for 
GRB21206. The IPN localized^ this GRB to 18° off solar, which precluded optical af- 
terglow searches; however, the brightness, duration, and proximity to the RHESSI rotation 
axis made it an ideal candidate to search for polarization. The shaded region shows our 5 s 
integration time for the polarization analysis. 
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Figure 2: The azimuthal scatter distribution for the RHESSI data, corrected for spacecraft 
rotation. Counts were binned in 15° angular bins between 0°-180°, and plotted here twice for 
clarity. The top plot shows the raw measured distribution (crosses), as well as the simulated 
distribution for an unpolarized source (diamonds) as modelled with a Monte Carlo code, 
given the time-dependent incident flux. The bottom plot shows the RHESSI data with the 
simulated distribution subtracted. This residual is inconsistent with an unpolarized source 
(dashed line) at a confidence level > 5.7a. The solid line is the best-fit modulation curve, 
corresponding to a linear polarization of (80 ± 20)%. 
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